DEVELOPMENT AND APPLICATION OF A TIME-HISTORY ANALYSIS 
FOR ROTORCRAFT DYNAMICS BASED ON A 
COMPONENT APPROACH 

Robert Sopher 

Supervisor, Dynamics Methods 

Daniel W. Hallock 
Dynamics Engineer 

Sikorsky Aircraft Division 
United Technologies Corporation, Stratford, CT 


Abstract 


This paper describes a time-history 
analysis for rotorcraft dynamics based on 
dynamical substructures, and non-struc- 
tural mathematical and aerodynamic com- 
ponents. The analysis is applied to 
predict helicopter ground resonance and 
response to rotor damage. Other applica- 
tions illustrate the stability and steady 
vibratory responses of stopped and gim- 
balled rotors, representative of new 
technology. Desirable attributes expected 
from modern codes are realized, although 
the analysis does not employ a complete 
set of techniques identified for advanced 
software. The analysis is able to handle 
a comprehensive set of steady state and 
stability problems with a small library of 
components. It has responded to new 
technologies with timely solutions bj 
limiting the effort required to implement 
new capabilities through its component 
structure . Opportunities were taken to 
reduce development costs by addressing 
more than one type of problem with a 
single component, such as using a minimum 
variance controller for trim and vibration 
reduction. 
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independent coordinate ^ 

system deunping matrix, 

N X N 

dependent coordinate [k]^ 

system damping matrix, 


rotor thrust coefficient 


D = matrix, Eq. (15), N x N 
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= damping matrix for i-th 
substructure 

= vector of external forces 
for i-th sub -structure, 
representing forces and 
moments . 

= components of force ap- 
plied to a substructure 
connection node, Fig. 2, 
lb 


= time interval in Newmark- 
Beta method, Eq. (10), 
sec. 

= independent coordinate 
system force Vector, N x 1 

= dependent coordinate 
system force vector, 

N^ x 1 

= matrix, Eq. (21), N x N 

= unit matrix, N x N 

= independent coordinate 
system stiffness matrix, 

N x N 

= dependent coordinate 
system stiffness matrix, 


= stiffness matrix for 
i-th substructure 

= independent coordinate 
system mass matrix, 

N x N 

= dependent coordinate 
system mass matrix, 


x mass matrix for the 
i-the substructure 
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components of moment 
applied to a substructure 
connection node, Fig. 2, 
lb- ft 

integer step number in 
time-history solution, 

Eg. (10) 

total number of indepen- 
dent coordinates in the 
assembled system. 

total number of coordi- 
nates in the dependent 
coordinate vector (ob- 
tained by summing the 
coordinates for all 
substructures ) 

matrix, Eg. (16), N x N 

matrix, Eq. (17), N x N 

vector of reaction loads 
applied to the i-th sub- 
structure 

components of reaction 
force applied to a 
substructure connection 
node. Fig. 2, lb 

components of reaction 
moment applied to a sub- 
structure connection node. 
Fig. 2, lb-ft 

matrix Eq. (13), N x N 

dependent coordinate 
system reaction load 
vector, Np x 1 

time, sec 

displacements at a con- 
nection node, ft 

vector of coordinates for 
i-th substructure 

rectangular coordinates, 
Fig. 2 

vector of dependent co- 
ordinates for the system, 
N 0 x 1, Eq. (3) 

vector of independent co- 
ordinates for the system, 

V x 1 

transformation matrix 
relating dependent coordi- 
nates to independent co- 
ordinates Eq. (4), Nx N 
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= Newmark-Beta factor 

= transformation matrices 
from which p is construc- 
ted, Eq. (9) 

- angular displacements at 
a connection node, 

Fig. 2, rad. 

= rotor solidity 


M 


= rotor advance ratio, non- 
dimensional 
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Subscripts 

n 


- rotor azimuth angle, rad 
= lotor speed, rad/s 
= matrix of zeroes 

= time index 


D = dependent coordinate 

system variable 

Superscripts 

i = i-th substructure 


X = transpose 

= first derivative with 
respect to time 

= second derivative with 
respect to time 


Introduction 


In the 1970's dissatisfaction with first 
generation computer programs for predict- 
ing helicopter performance and dynamic 
behavior motivated the development of the 
Second Generation Comprehensive Helicopter 
Analysis System (2GCHAS). The project is 
funded by the U.S. Government and is 
managed by the 2GCHAS Project Office at 
the NASA/Ames Research Center, and in- 
volves the participation of industry. The 
2GCHAS system aims to provide results for 
several helicopter related engineering 
disciplines, as well as helicopter dyna- 
mics. 

Several approaches were identified as 
being of potential value for overcoming 
first generation deficiencies. These 
approaches consisted of the use oi a 
unifying mathematical basis, executive- 
based software, and software design and 
management methodology. 
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One promising mathematical approach is to 
separate the dynamical structure into 
several components, or substructures, and 
subsequently to combine these into a 
system of second order differential 
equations. Coupled with automated 
assembly of che components, substructuring 
is expected to enable many problems to be 
modeled, overcoming the lack of 
versatility characteristic of first 
generation systems. Substructuring should 
reduce the difficulties of verifying the 
code by dividing .he system into easily 
verifiable parts. Activity can be focused 
on areas of new code during the process of 
adding new components, making the system 
more responsive to changes. 

In parallel with this, a software execu- 
tive would be used to enhance system 
versatility and usage for components which 
could not be handled as dynamical sub- 
structures, such as post-processing 
modules and certain types of aerodynamic 
components . 

Application of +ne techniques struc- 
tured design, developed in the software 
industry (Ref. 1), would help to improve 
architecture, and coding standards would 
make code legible. As a result, a second 
generation system would become more pro- 
ductive by being able to lead a 1 nnger 
useful life and would become more credible 
because it could be more easily verified. 
Finally, software management methodology, 
including automated software tools, would 
be uc »d for configuration control of 
versions of the system by protecting 
versions from untested and undocumented 
changes . 

The cost of applying all the approaches 
mentioned above to create a 3ystem with 
the scope of 2GCHAS appears to be beyond 
the resources of any single helicopter 
manufacturer. On the other hand, recog- 
nizing that benefits might result from 
application of a part of these approaches, 
the U.S. industry developed experimental 
codes utilizing some of the new concepts 
and limited to solution of dynamic pro- 
blems (Refs. 2 to 4). Sikorsky has 
developed two component-based codes con- 
sisting of the Simplified Vibration 
Analysis (SIMVIB) of Ref. 2, and the 
Rotorcraft System Dynamics Analysis 
(RDYNE ) described in this paper. Both 
methods utilize the same code for sub- 
structure assembly, but component 
libraries are different, providing 
solutions which differ in technique and 
scope. In contrast to SIMVIB which 
emphasizes harmonic balance solutionr, 
RDYNE is a time-history analysis. It i3 
the purpose of the paper to describe 
RDYNE, particularly with a view to its 
second generation attributes. 


Basis of Analysis 

The helicopter dynamical system is assumed 
to be made up of dynamical substructures 
and these are automatically assembled into 
a coupled system represent 2 d by a second 
order differential matrix equation. 
Figure 1 shows a typical substructure 
breakdown of a helicopter. The method is 
able to assemble components into many 
combinations and orientations . Coupled 
system response is obtained by integrating 
the differential equation with respect to 
time. 



Fig. 1 Sample Substructures Used 
in the RDYNE Rotorcraft 
Dynamics Analysis 


Substructure Assembly Method 

The substructure method employed is the 
Hurty method of Ref. 5. The coordinates 
or degrees of freedom of a substructure 
(also called physical component) are 
physical and generalised displacements, 
such as modal amplitudes. The matrix 
equation of motion fcr each substructure 
is expressed in mass, damping, and stiff- 
ness matrix, and force vector form. 

# ifhe a 

{r} ' . These are respectively mass, 

damping, and stiffness matrices, and 
external force and reaction force vectors. 
The submatrices for the substructures are 
collected into partitioned diagonal 
matrix equation which represents the 
system. This partitioned diagonal matrix 
equation is 

Vd + ®D*D + Vd = F D " *D (1> 
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The matrix Mq illustrates the typical form 
of the diagonal matrices as: 




tin] 


( 2 ) 



( 2 ) 


Substitution of the transformation Eq. (4) 
in Eq. (1) c'nd premultiplication of 
equation (1) by p 1 yields ti:e final equa- 
tion of motion for the coupled system. 
This is 

MXj + CXj + KXj = F (5) 


The vector of coordinates X- illustrates 
the form of the vectors F n and R-. in Eq. 
(1), and is d-d 

^ «/ { x} (1) 

(lxl (1 > 

The coordinate vector contains uncon- 
nected (or internal) ^coordinates and 
connection node coordinates . When com- 
ponents are to be connected to each other, 
redundant coordinates occur in X-. Fig. 2 
shows the translational and rotational 
displacements, and forces and moments on a 
connection node. 


The matrices in Eq. (5) are typified 

by T 

M = ^ i M D P (6) 

and the load vector is 

F = P T F d (7) 

The component displacements X.. may be 
recovered from the solution to°Eq. (5), 
X., by mapping Eq. (4), and component 
velocities and accelerations are derived 
similarly. With X D and its derivatives 
known, connection 0 node or interface 
reactions (e.g., rotor hub shears) may bt 
determined from Eq. (1). 

*D = F D - (M D*D + C A + W :8) 






Fig. 2 Displacements and Forces 
Acting on a Connection 
Node of a Substructure 


The synthesis of the equations of motion 
for the coupled system is accomplished by 
a napping relating the dependent co- 
ordinates, Xj., to a reduced (or indepen- 
dent) coordinate set, X,. Redundant co- 
ordinates are eliminated by requiring 
component displacements to be equal at 
connections. The transformation matrix 
relating X_ and X. is denoted by p, and 
the mapping of coordinates is 


The transformation matiix p is tla product 
of three transformations (Ref. 2) which 
assemble coupled systems from cmpon*. nts 
whose prope ties are defined in component 
local axes w.*d allow for the use of modal 
coordinates. The transformation p is 

P = P 0 Pl P 2 < 9) 


The transformation from dependent co- 
ordinates resolved to local axes with 
arbitrary angular orientations to depen- 
dent coordinates resolved to a global 
reference axis if p . The transformation 
from dependent coordinates, referred to a 
global reference axis, to coordinates from 
which redundant coordinates at cornection 
nodes have been removed is pj. The trans- 
formation from physical domain independent 
coordinates to coordinates which include 
modal coordinates, X x , is p 2 . 

Time History Solution 

The solution algorithm yielding the time 
history response is the Newmark-Beta 
method described in Ref. 6. Displacement 
responses are obtained from displacements 
known at prior tim^s , a*id from data 
defining the magnitudes of known external 
forces acting on the system. 

The vector of coordinates -atisfying Eq. 
(5) at time step n is denoted (X.) . The 
corresponding time is: 


*D * P X I 


(4) 
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In Eq. (10) h is the step size. The 

integer step number n ranges from one to 
the number of steps in the calculation. 
Time t is the initial time. Correspond- 
ing initial conditions are (X T r and 

<Vo' 

* 

ard (X.; are specified the 
responSe°at t^ is found from: 

(Xj)j = D' 1 [P(Xj) 0 + Q(X I ) o + 0*h 2 Fj + RFJ (11) 


When (X ) 
displacement 


Matrices .0, P, Q and R are functions of M, 
C, and K in Eq. (5), time step h, and 
factor 0* (ree equations (15) through 
(21)). They have dimensions NXN where N 
is the size of X^ in Eq. (5). The factor 
0* is the NewmarK-Beta Factor. Values of 
0 * ranging from 0 to 0.25 can be input to 
the program. The ractor 0* is used to 
control the variation of acceleration 
assumed in the time interval (see remarks 
in Ref. 6). A numerical solution with 0* 
=0.25 is unconditionally stable. Force 
vectors F 1 and F are known external 
forcing functions a? t Q and t^ . 

After the solution (X.). is known, suc- 
cessive solutions are 1 obtained from the 
recursion formula: 

«I>nn = D ' 1 lB 0 C I>n- F,( V„-l 


p 


M + (h/2)C - (1/2-0 *) h 2 K 




- ( 1/4-0 *)h3CM _1 K 

(16) 

Q 

= 

[M - (1/4 - 0*) h 2 CM -1 C] h 

(17) 

R 

= 

((1/2 - 0*)I + (1/4 -0*) 




hCff 1 ] h 2 

(18) 

I 

* 

unit matrix 

(19) 

B 

- 

2M - (1-20*/ h 2 K 

(20) 

F 1 

= 

M - (h/2 )C + 0*h 2 K 

(21) 

The 

dimensions of matrices in Eqs. 

(15) 

through (21) are NXN where N is 

the 


dimension of the vector of independent 
coordinates Xj. 

Main Features of Computer Program 

The top-level structure of the computer 
program reflects the component basis of 
the analysis. Dynamical components are 
substructures which obey second order 
differential equations which can be 
assembled into coupled dynamical systems. 
Non-dynamical components include several 
types of section aerodynamic and inflow 
models, and a component for trim and 
vibration reduction. 

Top-Level Structure of Program 


+ 0*h 2 (F q+1 + (l/p*-2)F n + F n _j)) (12) 

The NxN Matrices B and F* are known 
functions of M, C, and K, step size h, and 
factor 0* (see Eqs. (20) and (21)). 
Equation (12) may be used to restart 
solutions from solutions calculated up to 
the time of restart. 


The additional assumption is made in the 
program code that the force F - linearly 
extrapolates the forces F ana F , . For 
equal intervals of time h n this assumption 
reduces to: 


F . — 2 F — F -» 
n+1 n n-l 

end Eq. (12) becomes: 

<* 1 W D " 1 WVn- F,(X I ) n-l + h 2 r n 1 


(13) 

(14) 


Eq. (13) allows for the insertion of 
aerodynamic forces which are functions of 
X, and X T , by assuming F to be a function 
or Xj and its derivatives at a prior time. 

Matrices in Eqs. (11), (12), and (14) are: 


The top level structure of the RDYNE 
computer program is shown in Fig. 3 and 
reflects the component basis of the 
program, showing a separation into com- 
ponent-dedicated and component-independent 
areas. This organization is responsible 
for several desirable attributes. 
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D = M + (h/2)C + 0*h 2 K 


(15) Fig. 3 
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Component routines are required for each 
dynamical component listed in Table 1 and 
a typical organization is shown in Fig. 3. 
Component input and processing blocks form 
mass, damping, and stiffness matrices and 
aerodynamic and gravitational loads at the 
component level . Component- independent 
code performs assembly of substructures 
into a coupled system, reduction to 
independent coordinates, and solution of 
the equations . An interpretive routine 
(not shown in Fig. 3) reads names of 
selected components and corresponding 
input, and is followed by processing which 
utilizes component element numbers and 
connection node numbers (Fig. l) arid 
component orientations to assemble com- 
ponents into a coupled system. 

Non-dynamical components include several 
types of section aerodynamic and rotor 
inflow models, and a component for trim 
and vibration reduction. The components 
resemble dynamical components by having 
dedicated input and processing routines 
but otherwise do not behave like dynamical 
substructures because they cannot be 
assembled into dynamical systems by the 
assembly method described previously. 

Usage and understanding of the program are 
facilitated by ability of the program to 
limit input data to components selected by 
the user for his particular problem. This 
contrasts with first generation systems 
which required an understanding of the 
input for the most comp rehens xve problem 
even when preparation of only a part of 
this input was required. Processing and 


storage requirements are limited to the 
components selected, allowing resources to 
be tailored to each problem and to be less 
than the resources for the most comprehen- 
sive problem. Verification, modification, 
and addition of components is confined to 
component routines, making the system more 
responsive to change requirements. User 
and programmer experiences and data on 
execution speed have confirmed the above 
attributes . 

Dynamical Substructures 

All dynamical substructures, except modal 
structures (see text below) employ co- 
ordinates which are unconnected (or 
internal) coordixiates, and connection 
coordinates which enable the substructure 
to be as semi led with other substructures 
through the equating of displacements at 
the connection. With the exception of 
modal structures, it is necessary for each 
substructire to include six equations of 
equilibria'll corresponding to the six 
displacemer.ts at the connection node shown 
in Fig. 2, to enable it to be coupled to 
any other substructure. 

A noteworthy difference between dynamical 
substructures (Table 1) is that the blade 
models contain matrices which are explicit 
functions of time while the fuselage and 
matrix structure in Table 1 have constant 
matrices. Explicit time dependence occurs 
from a resolution of blade hub loads to a 
non-rotating axis system (Fig. 2) to 
derive the connection node equilibrium 
equations. This allows the transformation 


Table 1 - Dynamical Components in the RDYNE Analysis 


Component 
Elastic Blade 


Articulated Blade 


Modal Structure 


Matrix Structure 


Prescribed Force 


Fixed Absorber 


Description 

Normal modes elastic blade with flatwise, 
edgewise, and torsion elastic modes (Ref. 8), 
augmented to include six hub displacements, 
and expressed in M, C, K, and F forms. 

Simplified model containing a subset of the 
coordinates applicable to the elastic blade. 

Structure expressed in terms of normal 
mode coordinates. 

Generalized structure with fully populated 
M, C, and K matrices. 

Substructure providing for the application 
tc any component of a harmonically varying 
force of specified amplitude and phase. 

Vibration absorber which may be attached 
to any other substructure in the non-rotating 
system. 


158 


• i 




matrix p in Eq. (4) to be independent of 
time, and can be shown to justify the 
transformation of the dependent coordinate 
Eq. (1) to the independent coordinate Eq. 
( 5 ). 

The theory of Ref. 8, available in older 
codes at Sikorsky, was used as the basis 
of the elastic blade model (Fig. 4) to 
minimize the labor required to code and 
check *he model, and because this model 
has been substantiated through extensive 
comparisons with test data. Although the 
blade model is available in older codes , 
completely new coding was written to con- 
form to the component basis and legibility 
requirements of RDYNE. The equations for 
the internal coordinates were augmented 
with six hub equilibrium equations corre- 
sponding to hub connection node coordin- 
ates. The equations also were reduced to 
M, C, K, and F forms. 



Fig. 4 Schematic of Elastic Blade 
Substructure 


Table 2 - Mathematical and Aerodynamic Components 


Component 

Ti.ua History 
Integration Method 


Trim Controller 


Environment Input 


Aerodynamic Model 
Type 1 

Aerodynamic Model 
Type 2 

Aerodynamic Model 
Type 3 


Aerodyramic Model 
Type 4 


Description 

Component used to integrate the equations 
of motion with respect to time which 
employs the Newmark-Beta finite difference 
method of Ref. 6. 

Minimum variance controller used for rotor 
trim and coupled system vibration 
reduction (Ref. 10). 

Component defining the properties of the 
atmosphere. 

Simplified formula-based section aero- 
dynamic model for conventional airfoils. 

Table look-up section aerodynamic model 
for conventional airfoils. 

Simplified formula-based section aero- 
dynamic model for circulation control 
airfoils (Ref. 7). 

Table look-up section aerodynamic model 
for circulation control airfoils. 


Rotor Inflow Type l Momentum-based uniform rotor inflow 

component . 

Rotor Inflow Type 2 Variable rotor induced inflow component, 

using a matrix of wake influence co- 
efficients, calculated by the method 
described in Ref. 9 and transmitted tc 
RDYNE. 


Rotor Inflow Type 3 Momentum-based annulus inflow for hover 

applications. 

Rotor Inflow Type 4 Glauert inflow consisting of steady and 

first harmonic azimuthal and linear radial 
variations of .inflow. 
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The modal structure is an exception to the 
coordinate classification for other 
substructures. This substructure employs 
the normal modes of a substructure as 
coordinates and has constant mass damping, 
and stiffness matrices, and can be used to 
represent fuselages, and other systems, 
described in the applications. Each such 
structure is allowed to have up to five 
connection nodes, at which elements of the 
modal matrix for each r >rmal mode are 
defined, comprising three translations and 
three rotations, with the directions shown 
in Fig. 2. The substructure is coupled to 
others by expressing physical displace- 
ments at a connection as a summation of 
modes, to define P 2 , in Eg. (9). 


The aerodynamic generalized forces in the 
blade components are obtained from appli- 
cation of blade element theory and by 
invoking the section aerodynamic and 
inflow components listed in Table 2 to 
define aerodynamic properties. 

Aerodynamic Components 

Four section aerodynamic and four inflow 
models are available, and are listed in 
Table 2 . The section aerodynamic models 
include formula-based and table look-up 
methods. Input to the formula-based 
section Aerodynamic Model Type 1 consists 
of lift curve slope, maximum lift coeffi- 
cient, and coefficients required in 
expressions for drag and moment curves. 
This simple model is used when data on 
section characteristics are not accurately 
known, such as in blade damage simula- 
tions. Aerodynamic Model Type 2 is used 
to provide bivariant tables of charac- 
teristics expressed as functions of 
angle-of-attack and Mach number, obtained 
from wind tunnel tests. 

To link an aerodynamic component to a 
blade component, the user specifies in the 
input to the blade component the element 
number of the aerodynamic component, which 
is followed by corresponding input data or 
file names defining the location of data. 
The linkage procedure provides consider- 
able latitude for using different aero- 
dynamic components on different blades and 
blade sections, and has been well received 
by users. 

Rotor induced variable inflow is embodied 
in geometric influence coefficients trans- 
mitted to a file from an existing program 
external to RDYNE (Ref. 9). The wake form 
is a skewed helix and its geometry is 
assumed to depend on advance ratio. RDYNE 
determines by an iterative method wake 
circulations which are consistent with the 
inflow influencing blade section aero- 
dynamic loads. A strong coupling is 


achieved of variable inflow to rotor re- 
sponse in RDYNE without having to rely on 
a procedure involving iterative coupling 
of large programs, such as that used for 
SIMVIB in Ref. 2. 

Trim and Vibration Reduction 

A departure was made with older deriva- 
tive-based trim determination methodology 
by treating the trim problem as an opti- 
mization problem, and the opportunity was 
taken to solve both the trim problem and 
the vibration reduction problem. The code 
for the controller was derived from the 
vibration reduction program described in 
Ref. 10. 

A minimum variance controller was imple- 
mented which utilized an unconstrained 
minimization formulation to reduce dif- 
ferences between a target trim state, 
corresponding to specified steady hub 
loads, and components of steady hub loads 
for the actual rotor state. The transfer 
matrix relating hub loads to control 
inputs is initially calculated by a 
difference method from the results of 
perturbations to the control vector and 
subsequently is identified by a Kalman 
filter procedure, which is able to 
speedily identify the transfer matrix. A 
scalar performance index embodies the 
objectives to be minimized and differenti- 
ation of the performance index with re- 
spect to the control variables yields an 
optimal control state vector. Controls 
are updated according to the optimal 
formula at user-specified intervals. 

Figures 5 and 6 show the effects of 
application of the method to a four-bladed 
conventional rotor on a rigid support. 
The simulation employs a single elastic 
blade and multiplies the steady loa -v a 
factor of four to define the rotor ds. 
Figures 5 and 6 show that severe 'im 
objectives may be specified and .il- 
taneously achieved. 
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Fig. 5 Effects of Minimum Vari- 
ance Controller on Steady 
Hub LOw.ds For A Conven- 
tional Rotor 
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Fig. 6 Effects of Minimum Vari- 
ance Controller on Trim 
Control variables For A 
Conventional Rotor 


The same approach lends itself to vibra- 
tion reduction, with the vector of con- 
trolled variables containing coefficients 
of vibratory loads or accelerations at the 
blade passage frequency or at multiples of 
this frequency. Figures 7 and 8 show the 
effects of blowing controls on vibratory 
hub loads for a circulation control rotor, 
and the corresponding control state for 
selected blowing harmonics. A single 
blade analysis was used for this appli- 
cation although RDYNE is not limited to a 
single blade and can be used to reduce 
elastic airframe vibrations for a multi- 
blade rotor system. 



HUB 

LOAD 

LB 



Fig. 7 Effects of Minimum Vari- 
ance Controller on 4 Per 
Rev Hub Loads For Cir- 
culation Control Rotor 
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Fig. 8 Effects of Minimum Vari- 
ance Controller on Control 
Variables For Circulation 
Control Rotor 


A 



Fig. 9 Simultaneous Trim and 
Vibration Reduction on a 
Conventional Rotor 


The analysis is configured to nermit 
selection of rotor load or vibration 
objectives, singly on in different com- 
binations, and Figs. 9 and 10 show the 
combined effects of application of 
mechanical controls to achieve trim and 
vibratory hub load reduction. 



Fig. 10 Control Input History for 
Simultaneous Trim and 
Vibration Reduction on a 
Conventional Rotor 



It is seen that a comprehensive solution 
such as the use of a single minimum 
variance controller for trim and vibration 
reduction can be used to reduce develop- 
ment costs by addressing more than one 
type of problem. 


Applications 

Applications are described which show the 
comprehensive ability of RDYNE to handle 
stability and steady state vibratory load 
problems using a small number of elements 
and a time-history basis. Applications to 
new technologies are described. 

Ground Resonance Stability 

In early 1979 RDYNE was evaluated for its 
abiijty to predict ground resonance sta- 
bility involving the coupling of several 
articulated blades of the type shown in 
Table 1 to an elastic airframe. The 
elastic airframe was represented by a 
modal component (Table 1). The values of 
damping of the coupled systems were 
calculated from the decay in the responses 
of modal displacements in the airframe. 
FFT processing was not used. Figure 11 
compares the percentage damping predicted 
by RDYNE with results from a linear 
stability analysis (Ref. 11) and shows 
very satisfactory agreement. Figure 12 
compares stability boundaries from RDYNE 
with a Floquet solution (Ref. 12) for a 
rotor with a failed lag damper and again 
shows good agreement. The analysis is now 
used in a routine manner at Sikorsky to 
determine the stability of systems with 
failed lag dampers. It is seen that the 
time-history method can comprehensively 
treat problems which were solved by 
eigensolution and Floquet methods. 
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Fig. 11 Ground Resonance Stability 
Comparison Between RDYNE 
and Linear Stability 
Analysis 
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Fig. 12 Stability Boundaries 

Predicted by RDYNE And A 
Floquet Analysis Foi A 
Rotor With A Failed Lag 
Damper 


Non-symmetric Rotor and Damage Simulation 

The advent of RDYNE created new oppor- 
tunities for studying non-symmetric rotor 
and failure/damage simulations. 

The program was applied to a BLACK HAWK 
rotor to simulate the response of a 
rotor/drive train system to lightning 
induced damage to one of the blades. The 
effects of lightning were determined by 
passing electric currents through a 
portion of blade in ground tests. The 
aerodynamic characteristics of the blade 
after damage were estimated from its 
appearance and were loaded in as input to 
the simplified section aerodynamic model 
listed in Table 2. The normal modes of 
the drive train were inserted in the input 
to the modal component to represent the 
seven component drive train. Results from 
the program without blade damage were used 
as initial conditions for a restart 
solution with the damaged blade. Figures 
13 and 14 show blade and drive train re- 
sponses before and after the damage, A 
drive train schematic is also shown in 
Figure 14. The transition from four-per- 
rev to once-per-rev responses is evident. 
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MAIN ROTOR HUB 



TAIL ROTOR 



MAIN TRANSMISSION 


ANQLE GEAR BC <ES 


POWER TURBINES 



Fig* 13 Effects of Blade Damage on 
Time Histories of Blade 
Responses For A Coupled 
Rotor/Drive Train System 


Fig. 14 Drive Train Responses to 
Rotor Damage For A Coupled 
Rotor/Drive Train System 


The use of the modal component for the 
drive train in this problem and its use as 
an airframe in the ground resonance 
problems shows that versatility and 
reduced development cost can be achieved 
through using a single type of component 
for different types of problems. 

Stopped and Gimballed Rotors 

Stopped and gimballed rotors are new 
technologies to which RDYNE has been 
applied and for which it has provided 
timely solutions. 

Rotors which are stopped have no centri- 
fugal stiffening and have to be checked 
for bending divergence when in the swept- 
forward position. Figure 15 shows the 
increase with forward speed in deunping of 


the first elastic flatwise mode of a 
hingeleso rotor blade stopped at 135 
degrees azimuth. At the divergence speed 
the frequency of the mode is zero and the 
modal damping becomes infinite. The 
divergence speed agreed well with results 
predicted by the doublet lattice theory in 
the NASTRAN analysis. The ability of 
RDYNE to successfully model fixed lifting 
surfaces by means of an elastic blade 
component, originally derived for rotary 
wing applications, is the result of 
designing this blade component for mul- 
tiple applications. In contrast to the 
original derivation in Ref. 8, which 
non-dimensional ized variables by rotor 
speed, the new equations were left in 
dimensional form, allowing the non- 
rotating case to be treated without 
difficulty. 
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rig. 15 Divergence of a 45 Degree 
Forward- Swept Stopped 

Rotor Blade 



Fig. 16 General Arrangement of 
Dynaflex Rotor Hub and 
Schematic of Elastic 
Gimbal 


The Dynaflex gimballed rotor is a new 
concept under development at Sikorsky 
which incorporates unique features (Ref. 
13). Utilizing composite materials, the 
central hub is gimbal -mounted relative to 
the shaft achieving a universal joint 
action with a spring restraint to the 
tilting motion (Fig. 16). The arrangement 
is less cumbersome ttian a mechanical 
universal joint, and greatly reduces 


Coriolis forces. Lower blade loads and 
reduced vibrations should result. The hub 
design is an exceptionally clean aero- 
dynamic form. 


Tests have been conducted with a 4.4 ft 
radius model. Comparisons shown below are 
for a stiff inplane rotor tested in a 
hover rig from 350 to 650 rpm, and for a 
soft inplane rotor tested at advance 
ratios between .15 and .45, and C /a 
between .06 and .11. 


RDYNE was applied to this rotor to evalu- 
ate its stability and vibratory load 
bel avior. The analysis was modified m a 
period of two weeks to include the matrix 
substructure listed in Table 1 to simulate 
the gimbal. The short development was 
evidence o' the responsiveness of a 
component -be method to new technology. 


The rctf r/support system used in model 
tests was represented by means of four 
elastic blades, a matrix substructure with 
two connection nodes representing the 
gimbal , and a normal mode structure 
representing the support (Fig. 17). The 
upper connection node of the matrix 
substructure was attached tc the elastic 
blades, and the lower connection at the 
pivot bearing was attached to the normal 
mode support. The properties input to the 
matrix substructure were established by 
deriving gimbal mass, damping, and stiff- 
ness elements by the Lagrangian method. 
Stability results were obtained by means 
of a moving block method, and a single 
blade simulation was used to calculate the 
trim state of the gimbal rotor in level 
flight. 



Fig. 17 



Represent 

Rotor 


the 


MODAL 

SUPPORT 


to 
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Figure 18 compares test and analysis 
predictions of gimbal system frequencies 
and dampings in hover. Figure 19 compares 
test and analysis results for bending 
moments in level flight. Edgewise mode 
dampings established from a gauge on the 
moving blade are compared in Fig. 20. The 
hover stability comparisons are quite good 
as are flatwise bending moment variations 
with forward speed. The edgewise bending 
moments and stability in forward flight 
are overpredicted. 


M- 




I 









□ 
















□ 






OVMl 



IJ 

— 

— 








} 





n 







rJ 

td 

nz 







rerr 



M 

Ml 

MVt. 



) 












Had 

| 


r 






























































c 

IMU 


k. _ 


« 

a 


» A 

i 

1 








nrw 

i ? 

w — 

> < 

n 

■"H 

) < 


> < 

r~^ 

> < 







RMtUAMWOlU 

1 1 

5 



P 

H 

pa 

i ; 

i 




| 



N 

MW— V 
HAD 









i — 

L_ 


i 



L 




L- 

in 

5 







WOTOW — 


ftOVWV AMALTWS 
o TOT 



* ' I } } I T** * ' I I "I ‘I i 

w « 



Fig. 19 Comparison of vibratory 
Bending Moments for a 
Dynaflex Model Rotor in 
Level Flight. 


Fig. 18 Comparison of Predicted 
Frequencies and Dampings 
with Test Data For a Dyna- 
flex Model Rotor in Hover. 
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Fig. 20 Comparison of B? ade 
Edgewise Stability for a 
Dynaflex Model Rotor in 
Level Flight. 


Concluding Remarks 

An analysis utilizing dynamical substruc- 
tures and non- structural components has 
confirmed several desirable attributes 
expected fr a modern codes. These attri- 
butes include the ability to handle a 
comprehensive set of problems with a small 
library of components, supported by the 
ability to treat steady state vibratory 
and stability problems within a time- 
history framework. The analysis has 
responded to new technologies with timely 
solutions for advanced rotor concepts, by 
limiting the ef^rt required to implement 
capabilities through its component struc- 
ture. Comprehensive solutions, such as 
the use of a single minimum variance con- 
troller for trim and vibration reduction, 
have reduced development costs by address- 
ing more than one type of problem. 
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DISCUSSION 
Paper No. 11 



DEVELOPMENT AND APPLICATION OF A TIME- HISTORY ANALYSIS FOR ROTORCRAFT DYNAMICS 
BASED ON A COMPONENT APPROACH 
Robert Sopher 
and 

4 Daniel W. Hallock 



Marty Schroeder, Solar Energy Institute : You have talked about the modularity of your code and 

the flexibility of it — can you give some indication of the size of the code — RAM and ROM? 

Sopher : It’s actually 4 megabytes now. We are upgrading the IBM System to 3 megabytes so we 

are not that concerned about the size. We are not interested in overlay. At some point in time 
we may be interested in an executive which brings in only the routines that a user is interested 
in applying wMch would compile and link-edit the routines in object time to create a module 
which he was interested in using. So that's something we are interested in. 

Ed Austin, U.S. Army Applied Technology Laboratory : I have Just re^m^ of questions I would like 

to ask. I'll try to cut it down to Just a few. First, with regard to the aerodynamics. Have 
you given any consideration to the way your executive might handle aerodynamics that are not 
Just a function of the current state but are dependent on previous events, maybe values of the 
state vector or other dependent parameters? 

Sopher : We have given some consideration to that. For example, if you take the case of a 
general response in two-dimensional linear flow to an arbitrary impulse *ngle of attack change 
the response is provided by the Wagner function. In order to calculate the resulting lift you 
need to retain a history of what has happened to the motion and then you apply the kernel func- 
tion to that. It's a very simple thing, actually. So as far as I can see in that particular 
application all you h?' f e to do is store the history of the motion somewhere in the program. In 
regard to other types of application I haven't ally thought about anything other than that. 

Are you concerned about, say, a time-history representation of variable inflow? 

Austin : Something like that, yes. 

Sopher : [I have] not really thought about that very much. 

Austin : Another question regarding your algorithms for your controller. You showed all your 

variables as nice continuous functions of time. Do you actually treat them that way or do you 
only look at them once per revolution? 

Sopher : No. They go like this. They are discontinuous so really what I should have drawn was 

a set of . . . it would look something like a bar chart, but I Just drew a smooth curve through 
that. 

Austin : But what kind of algorithm is it that actually performs your convergence? 

Sopher : Well, the objective function is a weighted square of the things you want to minimize. 

For example, say, you want to achieve a certain level of lift. It is the square of the differ- 
ence between the target value of the lift and the actual lift as harmonically analyzed in the 
program, say, a steady value of lift. Now that is weighted by a weighting function. In addi- 
tion to that we have weights that are applied to the control vector as well and the purpose of 
that is to try to limit the magnitude of the excursions of the controls because this is an 
unconstrained minimum optimization approach. We haven't tried to use a constraint optimization 
approach because if you go to any of those programs like C0PES/C0 M MIN or ADS it's an incredibly 
large code and itself would equal the size of this program. So tt s was the simple approach 
[that was] adopted. The relationship between the change in the state of the variables which are 
being used to control the system and the controlled state is obtained through what we call the 
T-matrix. That’s identified subsequently by a Kaman filter method. As a matter of fact, I 
should make some acknowledgements here. Originally the controller was developed by John 
Molusis; Bob Taylor, who is with Boeing Vertol, went into further developments on it; an., 
then our research labs carried out further developments. We had to reconfigure it considerably 
before it could be usefully used in RDYNE. 

Austin : How often is that controller updated? 

Sopher : Well, as often as you like. The user has the ability to define how often he wants to 

do it, but I believe in applications that we have typically it's after each revolution or after 
2 revs or something like that. There are people who use the program who are more experienced in 
answering tnat question. It's undergoing fairly intensive use now. 

Austin : You show some Floquet results. Do you use the one-pass or the N-pass approach? 
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Sopher : No, those results were obtained by means of a time-history response. But I compared 

with Gene Hammond’s Floquet solution. We just use log decrement on that. In order to get the 
stability of the Dynaflex System we used a moving block method. As a matter of fact we haven't 
built the moving block into the program. Bob Goodman developed his own little post-processor. 

I think what would happen is you can eliminate some of the concerns about respoisiveness in 
terms of time for providing stability results by building a post-processor into the program. 

That would address some of the concerns people have about the comparative efficiencies of aero- & 

elastic stability methods versus time-history methods. 1 

Austin : One final question. Do you have your aerodynamics and dynamics integrated into one 1 

program or they separate programs coupled by JCL? 

Sopher : They are integrated into one program. 





